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Abstract 


This  report  documents  a  comprehensive  experimental  and  theoretical  investigation  of 
the  deformation  behavior  of  high-purity,  polycrystalline  <a-titanium.  A  series  of  mono¬ 
tonic  uniaxial  compression  and  tension  tests  were  carried  out  at  room  temperature  under 
quasi-static  conditions  to  quantify  the  plastic  anisotropy  and  the  tension-compression 
asymmetry  of  this  material.  The  evolution  of  microstructure  and  texture  during  deforma¬ 
tion  was  studied  using  optical  image  microscopy  (OIM)  and  neutron-diffraction  techniques 
to  elucidate  the  role  of  deformation  twinning  and  its  effect  on  the  macroscopic  response. 
In  addition,  four-point  bending  tests  were  conducted.  As  a  result  of  the  anisotropy  and 
directionality  of  twinning,  qualitative  differences  between  the  response  of  upper  and  lower 
fibers  in  different  orientations  were  observed.  To  characterize  the  material’s  strain  rate 
sensitivity,  Split  Hopkinson  Pressure  Bar  tests  at  strain  rates  of  400  to  600  sec-1  were 
performed  along  the  axes  of  symmetry  of  the  material.  A  clear  increase  in  strength  with 
increasing  strain  rate  is  observed,  the  hardening  rate  remaining  practically  unchanged 
for  all  directions,  with  the  exception  of  the  rolling  direction.  The  dramatic  hardening 
rate  increase  in  the  rolling  direction  was  indicative  of  higher  propensity  for  twinning  with 
increasing  strain  rate.  Based  on  the  experimental  data,  it  can  be  concluded  that  the  mate¬ 
rial  has  a  very  complex  anisotropic  behavior  and  exhibits  tension/compression  asymmetry 
and  strain  rate  sensitivity. 

The  experimental  activities  were  complemented  with  multi-scale  model  development 
in  the  framework  of  plasticity  theory.  To  describe  the  quasi-static  macroscopic  response  of 
the  material,  a  new  anisotropic  elastic/plastic  model  was  developed.  Key  in  its  formulation 
is  a  yield  criterion  that  captures  strength  differential  effects.  Anisotropy  was  introduced 
through  a  linear  transformation  on  the  Cauchy  stress  tensor  applied  to  the  material.  An 
anisotropic  hardening  rule  that  accounts  for  texture  evolution  associated  to  twinning  was 
developed.  It  was  demonstrated  that  the  proposed  model  describes  very  well  the  main 
features  of  the  quasi- static  response  of  high-purity  titanium  when  subjected  to  monotonic 
loading  conditions.  Validation  of  the  model  was  provided  through  comparison  between 
measured  and  simulated  strain  distributions  in  bending.  In  particular,  the  shift  of  the 
neutral  axis  was  well  described.  Rate  effects  on  the  behavior  of  the  material  were  also 
investigated.  This  integrated  experimental-theoretical  effort  appears  to  have  resulted  in 
the  description,  for  the  first  time,  of  the  anisotropic  stress-response  of  high-purity  a- 
titanium  at  room  temperature. 
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1  Introduction 


Because  of  their  outstanding  engineering  properties  such  as  corrosion  resistance,  200,  000  lb /in2 
strength,  moderate  weight  (40%  lighter  than  steel)  and  high  strengths,  good  formability,  and 
biocompatibility,  titanium  and  its  alloys  have  found  widespread  use  in  high-performance  appli¬ 
cations.  For  example,  titanium  and  its  alloys  are  widely  used  for  aircraft  structures  requiring 
greater  heat  resistance  than  aluminum  alloys  because  of  their  high  strength-to- weight  ratio. 
Also,  titanium  alloys  are  leading  candidates  for  metal-matrix  composites,  primarily  for  aircraft 
engine  applications.  The  stress-strain  response  of  titanium  at  high  strain-rates  and  its  alloys 
is  receiving  continued  attention,  which  is  related  to  platform  calculations  for  crash-worthiness 
and  foreign-object  damage  in  aerospace  systems,  ballistic  and  armor  applications,  high-rate 
forming,  and  high-rate  machining  [1].  The  defense  industry  is  also  interested  in  advancing  the 
application  of  titanium  alloys  to  armor  systems. 

Pure  titanium  has  a  hexagonal  close-packed  (HCP)  structure  with  c/a  ratio  of  1.587  lower 
than  the  ideal  c/a  ratio  of  1.633.  There  are  three  principal  types  of  titanium  alloys:  a  or  near-ct 
alloys,  a  —  (3  alloys,  and  (3  alloys.  Titanium  alloys  in  the  low  and  medium  temperature  regime 
mainly  consist  of  the  HCP  a-phasc  with  very  little  dispersed  (3  phase  in  between  the  a  grains. 
Considerable  efforts  have  been  devoted  to  the  characterization  and  modeling  the  response  of 
titanium  alloys  (see  for  example,  [2],  [3],  [4],  [5],  [6],  [7],  etc.).  The  experimental  data  reported 
concerned  mainly  the  uniaxial  and  torsional  responses  at  different  strain  rates  and  temperatures 
while  data  for  non-proportional  loading  paths  is  rather  limited  (see  [7]).  Many  studies  have 
also  been  devoted  to  the  understanding  of  the  specific  plastic  deformation  mechanisms  at  room 
temperature  in  pure  titanium  (e.g.  see  [1];  [8]  for  data  on  high- purity  a  titanium;  [9]  and  [10] 
for  commercial  purity  titanium).  Two  types  of  deformation  modes,  slip  and  twinning,  occur  in 
titanium  and  its  alloys  during  plastic  deformation  at  room  temperature.  It  is  generally  agreed 
that  pronounced  yield  asymmetry  is  associated  with  the  activation  of  twinning.  Previous  studies 
([7],  [11],  [13],  etc.)  have  also  shown  that  classic  plasticity  models,  such  as  J2  plasticity,  or  Hill 
[12]  are  unable  to  capture  this  asymmetry  that  results  from  the  combination  of  a  sharp  initial 
basal  texture  and  the  polarity  of  deformation  twinning.  Twinning  plays  two  important  roles  in  a 
titanium.  It  is  a  main  contributor  to  texture  evolution  by  reorienting  the  twinned  areas  (see  [10]) 
even  for  the  simplest  monotonic  loading  paths.  Furthermore,  twinning  drastically  influences 
the  strain  hardening  behavior  (see  for  example  [8]).  Accurate  and  realistic  modeling  of  the 
behavior  of  a  titanium  thus  requires  incorporation  of  the  effects  of  twinning  on  the  mechanical 
response.  Since  in  crystal  plasticity  models  the  distribution  of  crystal  orientations  in  the  given 
polycrystal,  the  available  slip/twinning  deformation  systems  and  the  stress  levels  necessary 
to  activate  them  are  taken  into  account  explicitly,  the  evolution  of  anisotropy  due  to  texture 
development  can  be  characterized  by  measuring  the  initial  texture  and  updating  the  texture 
using  a  suitable  homogenization  scheme  such  as  Taylor  or  self-consistent  model  (see  [14]). 
Recently,  the  application  of  crystal  plasticity  models  to  hexagonal  metals  and  the  incorporation 
of  crystal  plasticity  calculations  directly  into  finite  element  (FE)  analyzes  have  received  much 
attention.  Models  that  account  for  both  slip  and  twinning  activity  and  employ  Taylor  (e.g.  [15], 
[16]  or  self-consistent  ([17],  [18],  etc.)  averaging  schemes  to  predict  the  aggregate  behavior 
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have  been  proposed.  For  example,  [16]  concentrated  their  attention  on  the  development  of 
slip-hardening  and  twin-hardening  functions  based  on  the  experimental  data  on  high  purity 
titanium  deformed  in  simple  compression  reported  in  [8].  These  laws  were  incorporated  into  a 
Taylor  crystal  plasticity  framework  and  further  used  to  simulate  the  stress-strain  response  and 
texture  evolution  for  monotonic  loading  (uniaxial  compression  and  simple  shear).  However, 
no  attempt  was  made  to  predict  the  final  deformed  shape  of  the  specimens,  nor  to  perform 
benchmark  simulations  of  more  complex  monotonic  loadings  such  as  bending.  In  [19]  a  self- 
consistent  viscoplastic  model  linked  to  the  explicit  FE  code  EPIC  has  been  successfully  used 
for  describing  the  deformation  of  pure  zirconium  with  a  strong  initial  texture  under  quasi-static 
monotonic  loading  at  room  and  liquid  nitrogen  temperatures.  The  above  direct  implementations 
of  poly  crystal  models  into  FE  codes,  where  a  polycrystalline  aggregate  is  associated  with  each 
FE  integration  point,  have  the  advantage  that  they  follow  the  evolution  of  anisotropy  due  to 
texture  development.  However,  such  FE  calculations  are  computationally  very  intensive,  thus 
limiting  the  applicability  of  these  approaches  to  problems  that  do  not  require  a  fine  spatial 
resolution. 

An  alternative  promising  approach  is  to  develop  anisotropic  formulations  at  a  macroscopic 
level  that  can  be  easily  implemented  in  FE  codes,  and  thus  can  be  applied  routinely  for  de¬ 
tailed  analyzes  of  complex  forming  processes.  Unlike  recent  progress  in  the  development  of 
mathematical  descriptions  of  anisotropic  yield  surfaces  for  materials  with  cubic  structure  (for 
a  review  of  such  formulations  see  [20],  [21]  [22],  etc.),  macroscopic  level  modeling  of  yield¬ 
ing  and  strain-hardening  of  HCP  materials  is  less  developed.  Some  of  the  rigorous  methods 
proposed  to  account  for  initial  plastic  anisotropy  or  to  describe  an  average  material  response 
over  a  certain  deformation  range  (see  [23])  can  be  extended  to  HCP  materials.  The  major 
difficulty  encountered  in  formulating  analytic  expressions  for  the  yield  functions  of  hep  met¬ 
als  is  related  to  the  description  of  the  tension-compression  asymmetry  associated  to  twinning. 
Recently,  yield  functions  in  the  three-dimensional  stress  space  that  describe  both  the  tension- 
compression  asymmetry  and  the  anisotropic  behavior  of  HCP  metals  have  been  developed.  To 
describe  yielding  asymmetry  that  results  either  from  twinning  or  from  non-Schmidt  effects  at 
single  crystal  level,  Cazacu  and  Barlat  [24]  have  proposed  an  isotropic  criterion  expressed  in 
terms  of  all  invariants  of  the  stress  deviator.  This  isotropic  yield  criterion  was  applied  to  the 
description  of  crystal  plasticity  simulation  results  of  [25]  and  that  of  [26]  for  randomly  oriented 
polycrystals  (for  more  details,  see  [27]).  This  isotropic  criterion  was  further  extended  such  as 
to  incorporate  anisotropy  using  the  generalized  invariants  approach  proposed  by  [21].  Com¬ 
parison  between  this  orthotropic  criterion  and  data  on  magnesium  and  its  alloys  by  [30]  and 
titanium  (data  after  [11])  show  that  this  anisotropic  model  accurately  describes  both  anisotropy 
and  tension-compression  asymmetry  in  yielding  of  these  materials.  In  [29]  another  isotropic 
pressure-insensitive  yield  criterion  that  accounts  for  yielding  asymmetry  between  tension  and 
compression  associated  with  deformation  twinning  was  developed.  This  isotropic  criterion  in¬ 
volves  all  principal  values  of  the  stress  deviator.  To  account  for  both  strength  differential  effects 
and  texture-induced  anisotropy,  an  extension  of  this  criterion  using  one  linear  transformation 
operator  L  operating  on  the  Cauchy  stress  deviator  was  proposed  by  [29].  The  orthotropic  yield 
criterion  [29]  was  shown  to  exhibit  accuracy  in  describing  the  yield  loci  of  magnesium  (see  [28] , 


5 


[29])  and  titanium  alloys  ([29],  [6]  and  [7]).  A  macroscopic  model  that  captures  the  evolution 
of  anisotropy  due  to  evolving  texture  in  hexagonal  metals  subjected  to  monotonic  loading  con¬ 
ditions  was  proposed  by  [31]  .  Yielding  was  described  using  the  orthotropic  yield  criterion  [29]. 
To  model  the  change  in  shape  of  the  yield  locus  due  to  twinning  reorientation  of  the  lattice, 
evolution  laws  for  the  anisotropic  coefficients  involved  in  the  expression  of  the  yield  criterion 
[29]  were  developed  using  a  multi-scale  procedure  based  on  crystal  plasticity  calculations  of  flow 
stresses  with  the  self-consistent  viscoplastic  model  (VPSC)  model  [32]  and  macroscopic  scale 
interpolation  techniques.  Experimental  crystallographic  textures  and  stress-strain  curves  were 
used  as  input.  This  approach  was  successfully  applied  to  the  description  of  the  behavior  of 
zirconium  at  room  temperature  when  subjected  to  quasi-static  deformation  ([31])  or  dynamic 
loading  (  [33],  [34]) 

This  study  presents  a  comprehensive  set  of  experiments  aimed  at  understanding  the  response 
of  high-purity  titanium  at  room  temperature  and  consequently  modeling  the  observed  response 
at  different  strain  rates.  A  series  of  monotonic  quasi-static  uniaxial  tension  and  compression 
tests  were  conducted  to  quantify  the  plastic  anisotropy  and  the  tension-compression  asymmetry 
of  this  material.  In  addition,  four-point  bending  tests  were  carried  out.  To  characterize  the 
material’s  strain  rate  sensitivity,  Split  Hopkinson  Pressure  Bar  tests  at  strain  rates  of  400  to 
600  sec-1  were  performed  along  the  axes  of  symmetry  of  the  material.  To  gain  insight  into  the 
deformation  mechanisms,  OIM  and  texture  evaluations  were  also  conducted  in  conjunction  with 
the  mechanical  tests.  The  experimental  activities  were  complemented  with  multi-scale  model 
development  in  the  framework  of  plasticity.  The  material  parameter  identification  procedure 
will  be  outlined.  The  algorithmic  aspects  related  to  the  implementation  of  the  models  developed 
in  the  FEt  code  EPIC  will  be  presented  in  detail.  The  ability  of  the  proposed  model  to  capture 
the  main  features  of  the  observed  behavior  will  be  examined  by  comparing  the  experimental  data 
with  simulation  results  in  terms  of  stress-strain  response  in  tension/compression  and  bending 
for  different  orientations  with  respect  to  the  axes  of  symmetry  of  the  material. 


We  shall  use  the  following  notations.  The  inner  product  of  two  vectors  v  and  u  will  be 
denoted  by  v  ■  u.  The  dyadic  product  of  two  vectors  v  and  u  will  be  denoted  by  v  0  u. 
Superscript  T  indicates  the  transpose  operation  (A T)ij  =  Aji  for  any  second-order  tensor  A 
while  symbol  ltrl  denotes  the  trace  operator  applied  to  second-order  tensors,  i.e.  tr  (A)  =  Avl . 
The  symmetric  part  of  any  second-order  tensor  A  is  denoted  by  A.symm  and  is  defined  as 


Asymm  _  _  (a  _[_  a7)  ;  the  antisymmetric  part  of  any  second-order  tensor  A  is  denoted  by 


A.skew  and  is  defined  as  Askew  =  -  (A  —  AT) .  The  second-order  identity  tensor  is  denoted  by 

2 

I.  The  inner  product  of  two  second-order  tensors  A  and  B  is  defined  by  A  •  B  =  tr  (ATB), 
the  product  between  a  fourth-order  tensor  P  and  a  second-order  tensor  C  is  the  second-order 
tensor  defined  as  (C  :  P)y-  =  CijkiPki- 


6 


2  Experimental  procedures 


High-purity  o-titanium  was  mechanically  tested  at  room  temperature  under  both  quasi-static 
and  high  loading  rates.  To  gain  insight  into  the  deformation  mechanisms,  texture  measurements 
were  carried  out  in  conjunction  with  the  mechanical  tests  at  Los  Alamos  National  Laboratory 
(LANL).  Both  OIM  and  neutron  diffraction  techniques  were  used  to  evaluate  texture  evolution 
along  each  loading  path. 


2.1  Material 

The  material  used  in  this  work  was  high  purity  (99.999%)  titanium  that  was  purchased  from 
Alpha  Aesar  of  Johnson  Matthey  Electronics,  Inc.,  (Spokane,  WA).  The  material  was  supplied 
in  the  form  of  a  15.87  mm  thick  cross-rolled  disk  of  254  mm  diameter.  Optical  microscopy 
showed  that  the  as  received  material  has  equiaxed  grains  with  an  average  grain  size  of  about 
20  /rm  (Fig.  1). 


Figure  1:  Micrograph  of  the  high-purity  titanium  plate  material 

Twenty  samples  were  removed  using  water  jet  from  the  perimeter  of  the  plate.  Each  sample  is 
approximatively  19.05  mm  x  19.05  mm  x  19.05  mm  and  11.  32°  from  the  neighboring  sample  (see 
Fig.  2).  The  basal  plane  (0002)  figures  were  analyzed  to  determine  the  rolling  direction  (RD), 
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Rolling  Direction 


Figure  2:  Measured  pole  figure  of  the  basal  plane  (0002);  the  rolling  direction  (RD)  is  indicated 
in  red. 

which  was  subsequently  marked  on  the  plate  as  shown  in  Figure  2.  Note  that  the  plate  exhibits 
fiber  texture  resulting  from  the  rolling  process.  The  basal  plane  aligns  at  30°  to  40°  from  the 
plate  normal.  This  result  was  also  confirmed  by  the  analysis  of  the  [0001]  pole  figures  obtained 
from  neutron  diffraction  measurements  of  the  initial  texture  (see  Figure  3  )  carried  out  on  the 
neutron  time-of- flight  (TOF)  diffractometer  HIPPO  (High-Pressure-Preferred  Orientation)  at 
LANSCE  (Los  Alamos  Neutron  Science  Center). 

2.2  Quasi-static  behavior  of  high-purity  cr-titanium 

The  quasi-static  characterization  tests  consisted  of  uniaxial  tension  and  compression  tests  that 
were  conducted  at  room  temperature  and  at  a  nominal  strain  rate  of  0.001  sec-1,  using  an 
Instron  1125  testing  machine  .  Test  specimens  were  cut  from  the  plate  using  Electrical  Dis¬ 
charge  Machining  (EDM)  in  order  to  eliminate  artificial  hardening  that  may  occur  using  other 
machining  techniques.  Tensile  testing  was  carried  out  with  a  5  KN  load  cell  and  an  Instron 
model  number  G-51-12-A  extensometer  with  a  gauge  length  of  25.4  mm. To  characterize  the 
anisotropy  of  the  material,  standard  tensile  specimens  were  machined  such  the  tensile  direction 
was  either  parallel  to  the  rolling  direction  (RD),  transverse  direction  (TD),  or  45°  to  the  rolling 
direction,  respectively  (see  Figure  4  for  the  dimensions  of  the  samples). 
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ABC 


Figure  3:  (0001)  pole  figures  measured  by  neutron  diffraction  showing  the  initial  texture  in 
a-titanium  used  in  this  study.  The  rolling  directions  (RD)  and  transverse  direction  (TD)  are 
in  the  plane  of  the  disk,  while  TT  designates  the  the  through  thickness  direction,  normal  to  the 
plate.  In  A),  TD  direction  is  in  the  center  and  TT  direction  is  side  to  side,  in  B),  TT  direction 
is  in  the  center  and  TD  direction  side  to  side,  while  in  C),  RD  is  in  the  center  and  TD  side  to 
side. 

A  specialized  miniature  test  specimen  was  used  for  the  through-thickness  (TT)  tensile  tests 
(Figure  5).  In  order  to  examine  the  microstructural  evolution,  tests  were  carried  out  to  approx- 
imatively  10%,  20%,  30%,  40%  strains  respectively  or  until  complete  failure  of  the  specimen 
occurred.  Deformed  samples  were  thereafter  examined  by  OIM.  The  uniaxial  tensile  test  re- 
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Figure  4:  Geometry  and  dimensions  of  the  specimens  used  for  in-plane  uniaxial  tensile  testing. 

suits  (Figure  6)  show  that  the  material  displays  orthotropic  behavior.  An  anisotropy  ratio  for 
initial  yield  stress  defined  by  the  ratio  of  the  yield  stress  in  the  through-thickness  direction 
(the  largest)  to  that  in  the  transverse  direction  (the  smallest  value)  is  1.27.  In  all  tension  tests, 
shear-type  fracture  was  observed.  In  contrast,  tensile  fracture  of  magnesium  alloy  AZ31B, 
which  also  has  HCP  crystal  structure  is  brittle  (see  data  by  [35]).  To  examine  the  effect  of 
loading  orientation  on  the  mechanical  response,  cylindrical  compression  specimens  (  7.62  mm 
x  7.62  mm)  were  machined  such  that  the  axes  of  the  cylinders  were  either  along  an  in-plane 
(IP)  or  through-thickness  (TT)  plate  directions  (see  Figure  4).  Compression  testing  was  con¬ 
ducted  using  an  Instron  100  kN  load  cell  and  an  Instron  extensometer  model  number  G-51-17-A 
with  a  gauge  length  of  12.7  mm.  The  uniaxial  compression  results  (Figure  7)  show  an  initial 
yield  stress  anisotropy  ratio  of  1.18,  which  is  smaller  than  in  tension.  This  observed  variation 
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Figure  5:  Geometry  and  dimensions  of  the  through-thickness  uniaxial  tensile  specimen. 


in  compressive  yield  stresses  is  consistent  with  previously  observed  results  for  polycrystalline 
anisotropic  HCP  materials  such  as  AZ31B  magnesium  (see  [35])  and  high-purity  zirconium  (see 
[19]).  The  larger  compressive  flow  stress  in  the  through-thickness  direction  as  compared  to 
the  in-plane  orientations  is  due  to  the  basal  texture  of  the  material.  Under  through-thickness 
compression,  the  load  is  applied  essentially  perpendicular  to  the  basal  plane  thus  plastic  flow 
is  activated  for  higher  stresses  than  for  in-plane  samples,  which  have  more  favorable  conditions 
for  activating  prism,  pyramidal,  and  basal  slip.  Tensile  deformation  is  slip-dominated  hence  the 
anisotropy  in  yield  stresses  is  stronger  than  twinning- dominated  deformation,  which  is  observed 
in  compression.  Comparison  between  compression  and  tension  stress-strain  response  along  the 
rolling  direction  is  presented  in  Figure  8.  Although,  initially  there  is  no  significant  difference 
in  yielding  behavior  (at  0.02  %  strain  offset,  the  yield  stress  is  175  MPa),  a  very  pronounced 
tension-compression  asymmetry  is  observed  after  about  10%  strain.  Note  the  especially  sharp 
difference  in  hardening  evolution.  While  in  tension,  the  material  hardens  gradually  until  plastic 
localization  (necking)  occurs  at  about  30%  strain,  in  compression  strain-hardening  is  strongly 
non-linear,  with  a  very  pronounced  increase  in  hardening  rate  observed  at  10%  strain. This 
striking  increase  in  hardening  may  be  correlated  with  the  onset  of  deformation  twinning.  This 
hypothesis  was  verified  by  performing  OIM  measurements  of  the  deformed  specimens. The  OIM 
map  of  the  specimen  deformed  to  10%  strain  reveals  that  many  grains  have  twinned  (twins 
appear  red  in  Figure  9).  The  twin  volume  fraction  was  estimated  to  be  17%.  Figure  10  shows 
an  OIM  map  corresponding  to  20%  strain,  which  indicates  a  high  volume  fraction  of  twinned 
grains,  about  40%. No  other  loading  path  produced  this  level  of  twinning  activity.  These  results 
are  consistent  with  previous  observations  reported  by  [8],  [36]  on  polycrystalline  a-titanium 
of  similar  purity.  Furthermore,  as  in  the  case  of  the  material  investigated  in  this  study,  the 
maximum  twin  volume  fraction  was  observed  in  simple  compression  to  20%  strain,  the  reported 
volume  fraction  being  of  45%.  A  comparison  of  tension  versus  compression  response  along  the 
transverse  direction  is  shown  in  Figure  11.  Notice  that  the  stress- strain  curve  in  compression 
along  the  transverse  direction  does  not  show  the  features  present  in  the  stress-strain  response 
in  the  rolling  direction  compression.  No  significant  change  in  strain-hardening  rate  is  observed, 
which  is  indicative  of  minimal  deformation  twinning.  Indeed,  post-test  analysis  with  OIM  re¬ 
vealed  that  there  is  little  twining  activity  in  compression  up  to  40  %  strain  (less  than  5%  volume 
fraction) . 
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Figure  6:  Uniaxial  tensile  tests  results  along  the  rolling  (RDT),  transverse  (TDT),  and  through¬ 
thickness  direction  (TT) 


Comparison  between  quasi-static  test  results  in  monotonic  uniaxial  compression  and  tension 
along  the  through-thickness  direction  are  shown  in  Figure  12.  A  strong  tension/compression 
asymmetry  in  initial  yielding  is  observed  whereas  the  asymmetry  in  hardening  evolution  is 
less  pronounced.  Since  compressive  loading  is  applied  essentially  perpendicular  to  the  basal 
plane,  in  principle  deformation  twinning  should  be  favored  over  non-basal  slip.  However,  OIM 
measurements  have  shown  very  little  twinning  activity.  To  further  elucidate  the  role  of  the  de¬ 
formation  mechanisms  on  the  macroscopic  response  for  each  loading  direction  and  evaluate  the 
evolution  of  anisotropy,  post-test  texture  measurements  were  conducted  using  neutron  diffrac¬ 
tion. 


2.2.1  Microstructure  evolution  during  uniaxial  compression  deformation 

The  (0001)  pole  figure  of  the  specimen  deformed  to  10%,  20%,  30%  and  40%  strains  in  uniaxial 
compression  along  the  rolling  direction  are  shown  Figures  13  to  14.  A  significant  texture 
evolution  is  observed.  Note  that  at  20%  strain  the  c-axis  of  most  of  the  grains  has  flipped 
by  90  °,  which  correlates  with  the  OIM  post-test  analysis.  Figure  15  shows  these  same  pole 
figures  at  the  appropriate  strain  levels  on  the  stress-strain  curve.  It  is  clear  that  the  point  on 
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Figure  7:  Uniaxial  compression  tests  results  along  the  rolling  (RDT),  transverse  (TDT),  and 
through-thickness  direction  (TTT). 


the  stress-strain  curve  at  which  there  is  an  abrupt  change  in  hardening  rate  corresponds  to  the 
same  strain  level  at  which  a  strong  texture  evolution  is  observed.  Figures  16  and  17  show  the 
(0001)  pole  figures  for  the  specimens  loaded  in  uniaxial  compression  in  the  transverse  direction 
to  strains  of  10%,  20%,  30%  and  40%,  respectively.  Note  the  fairly  strong  alignment  of  the 
c-axis  in  most  of  the  grains  with  the  through-thickness  (TT)  direction.  Owing  to  the  basal 
texture  of  the  material,  the  compressive  flow  stress  in  the  through-thickness  direction  is  larger 
than  in  the  in-plane  orientations.  This  is  also  revealed  by  the  uniaxial  loading  tests  which 
show  that  the  plate  is  stronger  in  the  transverse  direction  than  in  the  rolling  direction.  Note 
that  there  is  only  a  slight  texture  evolution  indicating  that  deformation  twinning  may  not  be 
a  dominant  deformation  mechanisms  for  this  loading  path. 
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Figure  8:  Comparison  between  the  compression  and  tension  uniaxial  response  in  the  rolling 
direction. 


Figure  9:  Orientation  Imaging  Microscopy  map  showing  the  evidence  of  twins  (in  red)  in  the 
specimen  deformed  to  10%  strain  in  simple  compression  along  the  rolling  direction.  The  twin 
volume  fraction  is  17%. 
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Figure  10:  Orientation  Imaging  Microscopy  map  showing  significant  twinning  activity  (45% 
volume  fraction)  in  the  a-titanium  sample  deformed  to  20%  strain  in  simple  compression  along 
the  rolling  direction. 


Plastic  Strain 

Figure  11:  Comparison  between  the  compression  and  tension  uniaxial  response  in  the  transverse 
direction. 
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Figure  12:  Comparison  between  the  compression  and  tension  uniaxial  response  in  the  through¬ 
thickness  direction. 


10%  20% 

Figure  13:  Measured  (0001)  pole  figures  for  specimens  loaded  in  compression  to  10  and  20%  in 
the  rolling  direction  (RD  in  center  and  the  TD  direction  from  side  to  side). 
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Figure  14:  Measured  (0001)  pole  figures  for  specimens  loaded  in  compression  to  30  and  40  % 
strain  in  rolling  direction  (RD  in  center  and  the  TD  direction  from  side  to  side). 


Figure  15:  Texture  evolution  for  compressive  loading  in  the  rolling  direction 
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Figure  16:  (0001)  pole  figure  for  specimens  loaded  in  compression  to  10%  and  20%  in  the 
transverse  direction,  respectively  (TD  in  center  and  the  TT  from  side  to  side). 
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Figure  17:  (0001)  pole  figure  for  specimens  loaded  in  compression  to  30  %  and  40  %  strain  in 
transverse  direction;  TD  in  center  and  the  TT  from  side  to  side 


17 


Figures  18  and  19  show  there  is  only  a  slight  texture  evolution  in  uniaxial  compression  in  the 
through-thickness  direction  indicating  only  minimal  twinning  for  this  loading  path. 


001 


min  max 

-0.07  4.89 


10% 


001 


min 

-0.05 


max 

4  93 


20% 


Figure  18:  (0001)  pole  figure  for  specimens  loaded  in  compression  in  through  thickness  direction 
to  10  %  and  20  %  strain,  respectively;  TT  in  center  and  the  TD  direction  from  side  to  side 


30%  40% 


Figure  19:  (0001)  pole  figure  for  specimens  loaded  in  compression  in  the  through-thickness 
direction  up  to  30%  and  40  %  strain,  respectively;  TT  in  center  and  the  TD  direction  from  side 
to  side 
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Figure  20  and  Figure  21  show  these  same  pole  figures  at  the  appropriate  strain  levels  on 
the  stress-strain  curve  for  TD  compression  and  TT  compression,  respectively.  It  is  clear  that 
twinning  is  not  a  dominant  mechanisms  for  these  loading  paths.  There  are  no  abrupt  changes 
in  hardening  rate.  This  correlates  with  the  texture  measurements  that  indicate  minimal  grain 
rotations  during  deformation.  As  outlined  in  the  introduction,  the  mechanical  response  and 
deformation  mechanisms  of  high-purity  a-  titanium  have  been  reported  only  for  monotonic  sim¬ 
ple  compression  in  the  throuugh-thickness  direction,  plane  strain  compression,  or  simple  shear. 
The  texture  measurements  reported  here  expand  the  existing  knowledge  base  by  providing  a 
thorough  characterization  of  the  anisotropy  of  the  mechanical  response  and  its  evolution  during 
monotonic  loading.  Post  test  analysis  using  OIM  has  shown  that  significant  twinning  activity 
occurs  only  in  rolling  direction  compression. 


Figure  20:  Texture  evolution 


for  compressive  loading  in  the  transverse  direction. 
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Figure  21:  Texture  evolution  for  compressive  loading  in  the  through  thickness  direction. 
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2.2.2  Four-Point  Beam  Bending  Test  Results 


In  order  to  provide  validation  data  for  the  modeling  effort,  four-point  bending  tests  were  con¬ 
ducted. The  specimens  were  beams  of  length  57.15  mm  with  a  square  cross  section  of  6.35  mm 
as  shown  in  Figure  22  (A).  The  uniaxial  and  compression  tests  have  revealed  that  the  tension- 
compression  asymmetry  of  the  material  is  highly  directional.  Thus,  to  further  characterize  this 
tension-compression  asymmetry  by  bending  testing,  it  is  necessary  to  carry  out  tests  for  at  least 
four  test  configurations  and  loading  scenarios.  To  this  end,  four  beam  samples  were  cut.  Two 
samples  were  machined  such  that  the  long  axis  was  along  the  RD  while  the  other  two  samples 
were  machined  with  the  long  axis  along  the  TD.  For  each  of  these  two  sets  of  specimens,  one 
beam  was  loaded  in  the  through-  thickness  direction  while  the  other  sample  was  loaded  in  the 
in-plane  direction  normal  to  the  beam  axis.  These  four  test  configurations  are  shown  in  Figure 
22  (B).  To  simplify  the  notations  when  describing  the  beam  test  results,  we  denote  by  (x,  y,  z ) 
the  reference  frame  associated  to  the  three  orthotropic  axes  with  x  =  RD,  y  =  TD,  and  z  = 
TT.  The  testing  jig  with  a  test  specimen  is  shown  in  Figure  23.  The  two  upper  pins  were  dis- 
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Figure  22:  Four-point  bending  test  specimens:  A)  Specimen  dimensions;  B)  Orientation  def¬ 
initions:  in  Case  1  and  Case  2  the  long  axis  of  the  beam  is  aligned  with  the  rolling  direction 
(x) while  in  Case  3  and  Case  4  the  long  axis  of  the  beam  is  aligned  with  the  transverse  direction 

em¬ 
placement  controlled  to  approximately  5.5  mm.  A  typical  load  pin  displacement  path  is  shown 
in  Figure  24.  Along  one  side  of  each  beam,  a  speckle  pattern  was  sprayed  and  digital  image 
correlation  (DIC)  (see  [42,  43])  was  used  to  determine  the  strain  field  after  deformation. The 
image  taken  had  88  pixels  along  the  short  direction  of  the  beam.  The  beam  dimension  in  that 
direction  is  6.35  mm.  Thus,  the  physical  distance  between  pixels  is  72  micron/pixel.  The  DIC 
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Figure  23:  Four  point  beam  test  jig:  A)  with  test  specimen;  B)dimensions  and  pin  positions. 


method  can  detect  displacements  of  0.01  pixel,  therefore  the  error  is  less  than  1  /m i.  A  typical 
undeformed  and  deformed  grid  are  shown  in  Figure  25.  Note  that  the  grid  and  subsequently 
measured  strain  field  do  not  cover  the  entire  speckle  field.  Also, the  deformed  specimens  were 
cut  at  the  midpoint  along  their  axis  to  examine  the  final  deformed  cross  section.  Measurements 
at  this  cross  section  were  taken  for  comparison  to  the  FE  simulations.  A  contour  plot  of  the 
experimental  axial  strain  field  for  each  of  the  loading  scenarios  (see  Figure  22)  are  shown  in  Fig¬ 
ures  26  to  29. The  axial  strain  is  defined  as  the  strain  in  the  long  axis  direction  of  the  specimen. 
For  Case  1  and  Case  2  loadings,  the  long  axis  is  along  the  rolling  direction  x  so  the  axial  strain 
component  is  ex  while  for  the  Case  3  and  Case  4,  the  long  axis  corresponds  to  the  transverse 
direction  y,  therefore  the  axial  strain  is  ey.  In  all  tests,  some  non-uniform  deformation  occurred 
in  the  direction  normal  to  the  plane  for  which  the  data  were  reported,  thus  introducing  a  slight 
error  in  the  computation  of  the  axial  strains  using  the  DIC  methodology. 
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Figure  24:  Typical  load  vs  displacement  curve  in  bending:  the 
pin  displacement  (in  mm),  the  right  vertical  axis  represents  the 
sec)  is  plotted  on  the  horizontal  axis. 


left  vertical  axis  indicates  the 
applied  load  (in  kN),  time  (in 


Figure  25:  Typical 
relation  (DIC  )  for 


undeformed  and  deformed  beam  grid  pattern  used  with  digital  image  cor- 
generating  the  experimental  strain  field. 
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Figure  26:  Experimental  axial  strain  (sx)  fields  for  Case  1  configuration:  long  axis  in  x=RD, 
loading  in  y— TD. 
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Figure  27:  Experimental  axial  strain  (e^)  fields  for  Case  2  configuration:  Long  axis  in  x=RD, 
loading  in  z=TT. 
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Figure  28:  Experimental  axial  strain  {ey)  fields  for  Case  3  configuration:  Long  axis  in  y— TD, 
loading  in  x=RD. 
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Figure  29:  Experimental  axial  strain  (ey)  fields  for  Case  4  configuration:  Long  axis  in  y=TD, 
loading  in  z=TT. 
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3  Proposed  model  for  description  of  the  elastic-plastic 
deformation  of  high-purity  o-titanium 


To  model  the  behavior  of  the  high  purity  a-  titanium  material  in  the  quasi-static  loading 
regime,  an  elastic-plastic  modeling  approach  is  adopted. To  describe  the  onset  of  yielding,  an 
anisotropic  yield  criterion  that  captures  both  strength  differential  effects  and  the  anisotropy  in 
yielding  will  be  developed.  A  parameter  identification  procedure  based  on  experimental  results 
will  be  also  outlined. The  integration  algorithm  for  the  proposed  model,  and  its  implementation 
in  the  explicit  FE  code  EPIC  follows. Comparison  between  model  predictions  and  the  data  in 
uniaxial  loading  is  given  in  Section  3.5.  Validation  of  the  model  through  its  application  to  the 
description  of  the  bending  response  is  presented  in  Section  3.6. 


3.1  Proposed  yield  criterion 

One  of  the  goals  of  this  research  is  to  advance  the  current  state-of-the-art  by  developing  user- 
friendly,  micro-structurally  based  and  numerically  robust  macroscopic  constitutive  models  that 
can  capture  with  accuracy  the  particularities  of  the  plastic  response  of  hexagonal  metals,  in  par¬ 
ticular  high  purity  a-titanium.  A  full  three-dimensional  anisotropic  yield  criterion  is  proposed. 
Key  in  this  development  is  the  use  of  the  isotropic  yield  function  developed  in  [24]  that  cap¬ 
tures  tension/compression  asymmetries.  First,  a  brief  overview  of  this  yield  criterion  is  given. 
After  reviewing  the  general  aspects  of  linear  transformation  operating  on  the  Cauchy  stress 
tensor, the  anisotropic  yield  function  is  developed. The  input  data  needed  for  the  calculation  of 
the  anisotropic  yield  function  coefficients  are  then  discussed. 


3.1.1  Isotropic  yield  function 

If  the  deformation  mechanism  of  plastic  deformation  is  sensitive  to  the  sign  of  the  stress  as  is 
the  case  of  deformation  twinning,  the  isotropic  yield  function  ought  to  be  represented  by  an  odd 
function  in  the  principal  values  of  the  stress  deviator. To  describe  the  asymmetry  in  yielding, 
due  to  twinning,  Cazacu  and  Barlat  [24]  proposed  an  isotropic  yield  criterion  of  the  form 

/  =  4  -  cJz  =  4  (1) 

where  Ty  is  the  yield  stress  in  pure  shear  and  c  is  a  material  parameter;  J2  and  J3  are  the  second 
and  third  invariants  of  the  stress  deviator,  respectively.  The  constant  c  can  be  expressed  in 
terms  of  the  yield  in  uniaxial  tension,  crT,  and  the  yield  in  uniaxial  compression,  oc ,  as 

c=3V3(4-4>  (2) 

2  (crT  +  ac) 

If  <7t  =  ac ,  then  c  =  0,  and  the  isotropic  criterion  (1)  reduces  to  the  von  Mises  yield  criterion. To 
ensure  convexity,  c  is  limited  to  the  following  range:  c  E  [— 3\/3/2,  3\/3/4] .  In  principal  stress 
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(3) 


space  and  for  plane  stress  conditions,  the  yield  locus  is  expressed  as 
^  {ai  ~  ala2  +  5  -  ^  [ 2(J1  +  2(J2  ~  3(al  + 

When  c  7^  0,  this  yield  locus  is  a  triangle  with  rounded  corners.  As  an  example,  in  Figure  30 
is  shown  the  yield  locus  (3)  corresponding  to  three  different  ratios  of  <ttI&c ■  When  (Jt/ctc  is 
one,  the  yield  locus  reduces  to  the  von  Mises  ellipse. 


Figure  30:  Plane  stress  yield  loci  (3)  for  various  values  of  the  yield  stress  ratio  gt/ctc- 

Figure  31  shows  a  comparison  between  the  predictions  according  to  the  the  yield  criterion 
(1)  and  results  of  polycrystalline  simulations  for  randomly  oriented  FCC  (solid  circles)  and 
BCC  (open  circles)  polycrystals  obtained  by  Hosford  and  Allen  [25]  using  a  generalization  of 
the  Taylor-Bishop  model  [49].  In  the  polycrystalline  model  [25]  it  was  assumed  that  the  only 
deformation  mechanism  is  twinning.  According  to  this  polycrystalline  model,  for  FCC  randomly 
oriented  polycrystals  the  ratio  otI&c  is  0.78  which  corresponds  to  a  value  of  c  equal  to  0.92 
(see  (2  )  while  for  a  randomly  oriented  BCC  polycrystal  this  ratio  is  of  1.28,  hence  the  value 
of  c  is  —0.92.  It  is  clear  that  the  analytic  criterion  (1)  describes  very  well  the  asymmetry  in 
yielding  resulting  from  activation  of  twinning.  Next,  in  order  to  describe  both  the  asymmetry 
in  yielding  due  to  twinning  and  the  anisotropy  of  rolled  sheets,  an  extension  to  orthotropy  of 
the  isotropic  criterion  given  by(l)is  developed. 

3.2  Proposed  orthotropic  criterion 

An  extension  to  orthotropy  of  the  isotropic  criterion  (1)  can  be  obtained  by  using  a  linear 
transformation  approach.  This  method  consists  in  replacing  in  the  expression  of  the  isotropic 
criterion  (1)  the  Cauchy  stress,  cr  by  E  =  Lcr,  where  L  is  a  4th  order  tensor.  Thus,  the 
anisotropic  criterion  is  of  the  form 
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Figure  31:  Comparison  between  plane  stress  yield  surfaces  (with  cr3  =  0)  according  to  Cazacu 
and  Barlat  [24]  and  polycrystalline  model  of  Hosford  and  Allen  [25]  for  randomly  oriented  BCC 
polycrstals  (c  =  0.92);  and  FCC  randomly  oriented  polycrystals  (c  =  -0.92). 


f°  = 


(4) 


In  the  above  equation  ry  denotes  the  pure  shear  yield  stress.  The  tensor  L  satisfies:  (a)  the 
symmetry  conditions:  Ll3ki  =  Ljm  =  L]tik  =  LkUj  (  AbM  =  1-..3  ),  (b)  the  requirement  of 
invariance  with  respect  to  the  symmetry  group  of  the  material,  and  (c)  Lik  +  L2/c  +  L:>j.  =  0  for 
k  =  1,2,  and  3,  which  ensures  that  X  is  traceless  and  consequently  yielding  is  independent  of 
hydrostatic  pressure. Relative  to  the  orthotropy  axes  (x,y,  z) ,  L  is  represented  by 
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where  a*,  i=1...6  are  constants.  In  the  (x,  y.z)  frame,  x  represents  the  rolling  direction,  y  the 
transverse  direction,  and  z  the  thickness  direction.  Thus,  in  the  (x,  y,  z)  frame  the  transformed 
stress  tensor,  S,  is  given  by  : 
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Txz 


s  = 


g[(a2  +  a3)ax  -  a3ay  -  a2az]  a±Txy 

a4 7~xy  3  [  “t  (ai  "t  a3 )Oy  a^Tyz 

a5rxz  aeryz  l[~a2ax  -  aicry  +  (ai  +  a2)az]  . 


Let  J2°  =  [|tr  (E2)]  and  J3°  =  |  tr  (E3). 

Note  that  J2  is  an  extension  to  orthotropy  of  J2,  the  second  invariant  of  the  stress  deviator, 
while  J3  is  an  orthotropic  extension  of  J3,  the  third  invariant  of  the  Cauchy  stress  deviator.  In 
the  symmetry  frame,  ( x ,  y,  z ),  these  generalized  invariants  are  expressed  in  terms  of  the  Cauchy 
stress  cr  as 


j2  —  5  ( a2  +  a3  +  a2a3)  &x 
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J3  =  57  [  (a?a3  +  a2af)  <r3  +  (a?a3  +  a3ai)  a3  +  (a?a2  +  a3a?)  <r3 

+  (— aia?  +  aia?  —  a?a3  —  2 a2a|)  cr2ay 

+  (  aia2  —  aia3  —  a?a3  —  2 a2a3)  cr?crz 

+  (— a?a2  —  a?a3  +  a2a2  —  2aia3)  uyax 

+  (  a?a2  —  2a?a3  —  aia3  —  a2a3)  ayaz 

+  (— a?a2  —  a?a3  —  2aia?  +  a?a3)  cr?<7x 
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+  [a3a§<7x  -  (aia?  +  a3af)  ay  +  aiajaz]  t2z 

+  [(-a2a§  -  «3«i)  crx  -  agagcr^  +  a2a?crz]  r2z  } 

+  2tt4tt5  OjQTxyTXz^~yz 


For  plane  stress,  the  expressions  of  these  generalized  invariants  reduce  to 


T° 

J2 


g  [  (a2  +  a3  +  a203)  cr2  +  (a2  +  a2  +  aia3)  a2 
(—2a?  +  aia2  —  aia3  —  a2a3)  axay  ]  + 
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(7) 
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J3  —  -fj  [  (®2a3  d~  a“2a^)  & x  T  +  fl^ttl)  Cr3  + 

+  (— cqa^  +  aia2  —  a2a3  —  2a2a3)  crjjcq, 

+  (— afa2  -  a2a3  +  a2a\  -  2axa 3)  a2yax  ] 

+  ^  (  a2al<jxT2y  +  a\a\ayTxy)  (8) 

Hence,  the  proposed  orthotropic  yield  condition  is 

( JJ)3/2  -  cJJ  =  4  (9) 


The  anisotropy  coefficients  a,  and  the  strength  differential  parameter  c  involved  in  the  proposed 
anisotropic  criterion  (9)  can  be  determined  based  on  unixial  tensile  and  compressive  tests  along 
different  orientations  6  with  respect  to  the  rolling  direction  x.  In  particular,  if  Rtens  and  RCOmp 
denote  the  yield  stress  in  tension  and  compression  along  the  rolling  direction,  respectively,  then 
according  to  (9) 


p  _  Z21 

1  Hens  o 


(g-2  +  +  ^2^3)  2  —  C  (^2^3  a3a2)  j" 


(10) 


R, 
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comp 


(11) 


Similarly,  if  Ttens  and  Tcomp  denote  the  yield  stress  in  tension  and  compression  along  the  trans¬ 
verse  direction,  y,  then  the  model  predicts  that 


T, 


tens 


—  ^  (^1  +  ^3  ^1^3) 2  —  c  (alas  +  &3&i) 

>}■ 


Tcomp  —  |  (al  +  a3  +  ala3)  2  +  C  (^1^3  +  a3al 


(12) 


(13) 


When  <7i  =  a2  =  aj:  and  <r3  =  0  (equibiaxial  loading),  according  to  the  proposed  model  yielding 
occurs  when 


T 

ab  =ty 


(2ai  -  2b2  -  63)3/2  -  c 


&X&2  (<2l  +  Q.2) 
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(14) 


rc 


while  under  equibiaxial  compression,  i.e.,  for  ax  =  a2  =  a'h 

CLid2(di  +  ci2) 


a?  =  r„ 


(2ai  -  262  -  63)3/2  +  c- 
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(15) 


The  anisotropy  coefficients  involved  in  the  proposed  orthotropic  yield  criterion  can  be  deter¬ 
mined  by  minimizing  cost  functions.  The  experimental  data  in  the  cost  functions  may  consist 
of  yiled  stresses  in  tension  and  compression  corresponding  to  different  orientations  in  the  plane 
of  the  plate  and  normal  to  the  plane  of  the  plate  using  the  above  equations.  If  available,  r- value 
data  can  be  used.  If  experimental  data  are  not  available  for  a  given  strain  path,  they  can  be 
substituted  with  numerical  data  obtained  from  polycrystalline  calculations. 
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For  the  o-titaniuin  material  studied,  the  coefficients  of  the  model  corresponding  to  fixed 
levels  of  accumulated  plastic  strain  (up  to  0.5)  were  determined  using  the  measured  tensile  flow 
data  at  0°,  45°,  90°  and  compressive  flow  data  at  0°  and  90°.  In  tensile  loading,  the  specimens 
experienced  uniform  deformation  up  to  about  20%  strain,  at  which  level  localized  deformation 
and  necking  occurred.  Tensile  data  beyond  0.2  strain  were  obtained  by  extrapolation.  The 
corresponding  theoretical  yield  surfaces  along  with  the  experimental  values  (filled  squares) 
are  shown  in  Figure  32.  Note  that  the  proposed  criterion  matches  the  data  very  well.  The 


Figure  32:  Theoretical  yield  surfaces  according  to  Equation  (9)  and  experimental  data  (symbols) 
corresponding  to  various  strain  levels 


quadratic  yield  criterion  of  Hill  [12]  is  the  most  widely  used  orthotropic  yield  criterion  available 
and  has  proven  to  be  accurate  and  robust  for  many  materials,  especially  steels.  However,  it 
can  not  account  for  the  strength  differential  observed  in  hexagonal  materials.  For  comparison 
purposes,  Hill’s  criterion  is  applied  to  the  high  purity  titanium  material  used  in  this  research. 
With  respect  to  the  orthotropy  axes  (x,y,z),  Hill  [12]  orthotropic  yield  criterion  is  written  as 

F(ay  —  az)2  +  G(az  —  crx )2  +  H(ax  —  ay )2  +  2  Lt2z  +  2  Mrzx  +  2  Nrxy  =  1  (16) 

where  the  coefficients  F,  G,  H,  L,  M,  and  N  are  material  constants.  The  quasi-static  ex¬ 
perimental  data  on  high-purity  a  titanium  plate  given  in  Table  1  and  Table  2  were  used  to 
determine  the  coefficients  (given  in  Table  3)  of  this  criterion. 
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Table  1:  Compressive  yield  data  for  high-purity  a-titanium  used  to  identify  Hill  [12]  parameter 

values. _ 

Direction  x  y  z 

Yield  Strength  (MPa)  142.7  208.5  246.8 


Table  2:  Tensile  yield  data  for  high-purity  a-titanium  used  to  identify  Hill’s  [12]coefficients. 
Direction  x  45°  y  z 

Yield  Strength  (MPa)  127.1  148.5  200.8  255.1 


In  Figure  33  are  shown  the  theoretical  yield  loci  according  to  Hill  [12]  along  with  the  yield 
loci  predicted  by  the  anisotropic  model  developed  as  part  of  this  research,  and  data  on  high- 
purity  a-titanium.  Note  that  Hill’s  [12]  criterion  cannot  capture  the  observed  behavior  while 
the  proposed  model  describes  very  well  the  pronounced  strength  differential  effects  displayed 
by  the  material. 


3.3  Modeling  anisotropic  hardening  during  monotonic  deformation 
of  ct-titanium 


From  the  experimental  data  (see  Figure  32,  it  appears  that  there  is  distortion  of  the  shape  of 
the  yield  surface  of  a-titanium  even  for  the  simplest  monotonic  loading  paths.  This  anisotropic 
hardening,  which  is  due  to  evolving  texture,  cannot  be  described  with  existing  isotropic  or 
kinematic  hardening  laws.  Recently,  [31]  proposed  a  method  for  accounting  for  texture  evolu¬ 
tion. This  method  allows  for  the  variation  of  the  anisotropy  coefficients  with  accumulated  plastic 
deformation,  i.e.  the  linear  transformation  operator  L  is  no  longer  constant. However, obtaining 
analytic  expressions  for  the  evolution  laws  of  all  the  Ly-  components  would  be  a  difficult  task. 
An  alternate  method  is  thus  proposed.  Specifically,  using  experimental  yield  stresses  and  nu¬ 
merical  test  results,  the  coefficients  of  the  anisotropy  tensor  L  involved  in  the  anisotropic  yield 
criterion  (9)  are  determined  for  a  finite  set  of  values  of  the  equivalent  plastic  strain,  e.  Next  are 
calculated,  the  effective  stress  defined  according  to  the  proposed  criterion  and  Y:1  =  Y(cr:! ,  Y), 
where  Y  defines  the  effective  stress-effective  plastic  strain  relationship  in  a  given  direction  (e.g., 
the  rolling  direction).  Further,  an  interpolation  procedure  is  used  to  obtain  the  yield  surfaces 
corresponding  to  any  given  level  of  accumulated  strain.  Indeed,  note  that  the  effective  stress 
according  to  the  proposed  criterion  can  be  written  as 


a  =  A\ 


(17) 


where  A\  is  a  constant  defined  such  as  to  assure  that  a  reduces  to  the  tensile  yield  stress  in  the 
rolling  direction. Hence, 
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Table  3:  Hill’s  [12]coefficients  for  high-purity  a-titanium. 
Hill’s  (1948)  Coefficients  F  G  H 

Value  2.34E-05  -7.598E-06  3.15E-05 


Figure  33:  Theoretical  yield  loci  according  with  Hill  (1948)  criterion  and  the  proposed  criterion 
9  in  comparison  with  experimental  flow  stresses  (symbols)  on  high-purity  a-titanium. 


A\  — 


-  ((^2  +  a3  +  a2®3)  ^  —  ^®2^3(a2  +  ^3) 


-1-1/3 


(18) 


For  an  arbitrary  value  of  the  accumulated  plastic  strain  <  £  <  £J+1,  the  anisotropic  yield 
surface  is: 

^current (*T  ^)  ^current  (  rr .  c  ) 

where 

d(cr,  §)curreni  =  £  <xJ  +  (1  -  £)crJ+1 

and 

the  interpolation  parameter  £,  being  defined  as: 


£  =  m 


ej+1  -  e 


£j+1  -  ej 


To  simulate  the  response  of  high-purity  a-titanium,  we  will  use  the  orthotropic  criterion  given 
by  (4)  and  associated  flow  rule  with  hardening  as  described  by  the  interpolative  procedure 
described  above. 
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3.4  Integration  algorithm  of  the  proposed  model,  finite-element  im¬ 
plementation,  and  simulation  of  the  response  of  high-purity  a 
titanium  in  uniaxial  tension/compression 

The  developed  orthotropic  yield  criterion  (see  (9))  in  conjunction  with  the  anisotropic  hard¬ 
ening  law  presented  in  the  previous  section  and  associated  flow  rule  was  implemented  into 
EPIC2003,  the  2003  version  of  the  explicit  finite  element  code  EPIC  (Elastic  Plastic  Impact 
Calculations). EPIC  code  has  been  developed  by  Dr.  Gordon  Johnson  under  the  primary  spon¬ 
sorship  of  the  U.S.  Air  Force  and  U.S.  Army.  The  first  documented  version  proposed  in  1977 
was  2D  only  but  has  later  evolved  into  a  one,  two  or  three-dimensional  version  with  many 
additions  and  enhancements.  All  simulations  were  carried  out  on  a  PC  platform  using  Compag 
Visual  Fortran  Professional  Edition  6. 6. a.  The  code  was  compiled  such  that  all  real  variables 
were  double  precision. 

Given  the  current  stress  state,  the  current  strain  rate,  and  the  current  time  step,  a  user’s  sub¬ 
routine  was  developed  in  order  to  update  the  stress  state. The  integration  algorithm  solves  for 
the  updated  stresses  by  enforcing  the  Kuhn-Tuckner  and  consistency  conditions.  As  already 
mentioned,  an  associated  flow  rule  was  assumed  such  that  the  stress  potential  G  is  the  same  as 
the  yield  function  f°  (see  (9)) 

Aep  =  A\J-  (19) 

where  Aef)  is  the  plastic  strain  increment,  A  A  is  the  plastic  scalar  multiplier,  f°  is  the  pro¬ 
posed  anisotropic  yield  function,  and  cr  is  the  Cauchy  stress  tensor. Within  the  stress  update 
subroutine,  an  updated  trial  stress  is  computed  assuming  elastic  response,  i.e.  dc r  =  Cedc, 
where  Ce  is  the  fourth-order  elastic  compliance  tensor.  The  updated  trial  stress  is  computed  as 
a  updated  =  &  current  +  da.  Next,  an  equivalent  stress  (a)  is  computed  according  to  the  proposed 
anisotropic  yield  criterion  and  compared  to  the  current  value  of  the  hardening  function  (I5). 
If  a  <  Ys  the  current  stress  state  is  elastic,  the  trial  stress  coincides  with  the  actual  stress.  If 

a  >  Ys,  then  an  implicit  iterative  scheme  is  used  to  return  the  stress  state  to  the  yield  surface. 

The  initial  guess  for  this  iterative  process  is  a.  The  total  accumulated  plastic  strain  (e)  is  taken 
as  the  hardening  parameter,  the  yield  criterion  being 

f(a,e)  =  a(a,e) -Ys(e)  (20) 

In  the  following  j  is  used  as  the  global  counter,  i.e.  j  is  the  current  state,  and  j  +  1  is  the  fully 
updated  state  and  n  is  used  as  a  local  iteration  counter,  i.e.  n  is  the  current  iteration  and  n  + 1 
is  the  updated  iteration  count.  Then,  the  total  increment  of  stress  for  a  given  iteration  step  is 

=  Aoi+1  +  S(TnX\  (21) 

where  A  indicates  an  increment  over  the  entire  time  step  and  5  indicates  an  increment  for  each 
iteration.  The  plastic  correction  to  the  stress  is 


The  derivatives  of  the  equivalent  stress  are  found  by  taking  only  the  first  term  of  a  Taylor 
series  expansion  about  the  current  state 
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Keeping  only  the  first  term  in  (23),  Equation  (22)  becomes 

x<Ti+1  ~  ni+!  'J+1 

0(7 n+ 1  ~  0An+ 1  L  ■ 

Solving  for  the  increment  for  a  given  iteration  gives 


X-.7+1  =  _r,~L  ■  — 

^^77+1  ^  *  ^^77+1^^. 


J  +  l 


(24) 


(25) 


A  Taylor  series  expansion  of  the  yield  criterion  is  used  to  obtain  an  approximation  of  the 
increment  of  the  effective  plastic  strain 


/+ 


j+ 1 

77+1  ?  ^77+1  / 


f(a^+1  -\ — — 

J  \u  n  )  '  q- 


3+ 1 


«;ii = o 
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Evaluating  derivatives  at  the  previous  step,  equations  (25)  and  (26)  can  be  used  to  give 


Wli  = 


/K+yy1) 

dd_  _i_  dYs  _  da_ 

da  °  da  ^  de  de 


(27) 


This  can  now  be  used  in  equation  (25)  to  find  8a^l  and  finally,  the  total  stress  increment 
is  found  from  equation  (21).  This  new  stress  is  then  evaluated  to  see  if  it  has  converged  within 
a  specified  tolerance.  If  not,  this  is  used  as  the  starting  stress  for  the  next  iteration.  When  the 
stress  has  converged,  the  global  stress  tensor  is  updated  and  returned  to  the  main  program. 
Derivatives  of  the  yield  function  with  respect  to  the  stress  are  used  in  integrating  the  increments 
of  plastic  strain  in  the  FE  implementation.  The  general  form  can  be  written  as 


dF  _  dF  dJ%  dS  dF  d  J3°  dS 
da  ~  dJZ&Zda  +  d.JS  dS  da 

or 

dF  _  dF  dJ%  0Zkl  dF  dJ%  dFkl 
daij  dJ°2  dT,kt  daij  +  dJ£  dT,kt  da^ 

where 

OF  _  0  (+)1/2 

dJ2  2  ((J2°)3/2  -  c  +)2/3 


35 


dF  _  0  c 

~  3  (( J2°)3/2  -  c  J|)2/3 

The  derivatives  with  respect  to  S  are 
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All  other  derivatives  are  equal  to  zero. 

The  FE  implementation  uses  the  methodology  for  describing  the  evolution  of  the  yield  surface 
due  to  texture  changes  presented  in  the  previous  section.  A  sufficient  number  of  points  was 
used  to  ensure  that  the  hardening  curve  was  recreated  with  enough  accuracy.  For  each  of  these 
strain  levels,  ranging  from  0  to  50%,  the  corresponding  yield  surfaces  according  to  Equation 
(9)  were  determined.  The  anisotropy  coefficients  shown  in  Table  4  were  calculated  following 
the  identification  procedure  described  in  Section  3.2. 

For  the  material  investigated,  the  iso-strain  contours  of  the  theoretical  biaxial  yield  surfaces 
compared  to  data  are  shown  in  Figure  34  A.  The  data  used  to  identify  the  model  parameters 
are  represented  in  Figure  34  B  by  symbols.  Recall  that  to  obtain  the  tensile  data  for  strains 
larger  than  20%  strain  required  an  extrapolation  because  the  material  began  to  localize  beyond 
this  strain  level.  Note  that  the  model  reproduces  the  data  quite  well  with  the  largest  difference 
between  theoretical  and  experimental  values  being  for  the  biaxial  stress  state.  Table  4  gives 
the  uniaxial  stresses  in  the  RD  and  anisotropy  coefficient  values  determined  at  each  strain  level 
for  our  material. 
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A  B 

Figure  34:  A)  Theoretical  yield  surfaces  for  fixed  levels  of  accumulated  plastic  strain;  B)  Data 
used  in  identifying  the  theoretical  yield  surfaces. 

3.5  Application  of  the  proposed  model  to  cr-titanium 

The  proposed  model  was  used  for  simulation  of  the  response  of  the  cc-titanium  material  inves¬ 
tigated  for  uniaxial  loading  condition.  The  model  was  then  validated  by  simulating  the  four 
point  bend  tests  described  in  the  section  3.6.  A  comparison  was  made  in  a  qualitative  way 
by  juxtaposing  contours  of  the  experimental  data  against  the  results  from  the  simulation  for 
all  four  loading  scenarios.  Comparison  of  the  experimental  cross  sections  of  the  beams  and 
simulated  ones  using  the  model  and  an  isotropic  material  with  a  von  Mises  yield  surface  were 
also  performed. A  quantitative  comparison  was  done  by  comparing  axial  strain  versus  height  at 
the  centerline  of  each  beam  specimen. 


Table  4:  Anisotropy  coefficients  for  discrete  strain  levels 


Strain 

Yield 

al 

a2 

a3 

a4 

a5 

a6 

c 

0.0000 

208 

0.5454 

0.5010 

1.0900 

0.7246 

-0.8675 

-0.8675 

-0.2168 

0.0250 

245 

0.5231 

0.4745 

0.9034 

0.7309 

0.7202 

0.7202 

-0.2198 

0.0500 

261 

0.6694 

0.5585 

1.1030 

0.9138 

0.9381 

0.9381 

-0.2291 

0.0750 

273 

0.6960 

0.5969 

1.1270 

0.9838 

0.9716 

0.9716 

-0.2607 

0.1000 

284 

0.5356 

0.4768 

0.8603 

0.7761 

0.7714 

0.7714 

-0.2754 

0.2000 

317 

0.0610 

0.0576 

0.0869 

0.0870 

0.0794 

0.0794 

-0.5908 

0.4000 

370 

0.0632 

0.0620 

0.0788 

0.0816 

0.0801 

0.0801 

-1.0330 

0.5000 

389 

0.9547 

0.9570 

1.2140 

1.1810 

1.1760 

1.1760 

-1.1480 

Note:  Yield  Strength  in  MPa 
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For  uniaxial  loading  conditions,  the  simulations  involved  a  single  element  with  eight  nodes 
with  a  single  integration  point  (see  Figure  35). The  cell  was  stretched  uniaxially  and  the  ob¬ 
tained  stress  versus  strain  results  were  compared  to  the  appropriate  experimental  data. For  each 
simulation,  four  nodes  on  one  face  of  the  element  were  restrained  and  the  four  nodes  on  the 
opposite  face  were  given  a  constant  velocity  in  either  the  tensile  or  compressive  direction. 

The  simulated  normal  stresses  in  the  appropriate  loading  direction  versus  the  effective  strain 


Figure  35:  Single  cell  computational  configuration 

are  compared  to  the  corresponding  experimental  data  (see  Figures  36  to  38).  Note  that  the 
model  accurately  reproduces  the  data  for  each  loading  condition.  The  largest  discrepancy  oc¬ 
curs  for  the  90°  tensile  data,  which  is  due  to  the  slight  discrepancy  between  experimental  and 
predicted  flow  stresses  in  this  direction  that  was  noted  previously. 


0  10  20  30  40  50 

Effective  Strain  (%) 


A  B 

Figure  36:  Comparison  between  finite  element  simulation  results  using  the  proposed  model  and 
experimental  data  corresponding  to  uniaxial  loading  in  the  rolling  direction:  A)  tension;  B) 
compression. 


A  B 

Figure  37:  Comparison  between  finite  element  simulation  results  using  the  proposed  model  and 
experimental  data  corresponding  to  uniaxial  loading  in  the  transverse  direction:  A)  tension  B) 
compression. 
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Figure  38:  Comparison  between  finite  element  simulation  results  using  the  proposed  model  and 
experimental  data  corresponding  to  uniaxial  loading  in  the  through-thickness  direction  :  A) 
tension  B)  compression. 


3.6  Model  validation  results 

Finite-element  simulations  of  beam  bending  tests  were  conducted. The  simulated  final  cross 
sections  for  the  four  configurations  are  shown  in  Figure  39.  As  expected,  when  the  hard 
direction  (TT)  is  perpendicular  to  the  loading  direction  (Cases  1  and  3)  the  cross  section 
remains  nearly  square.  Case  2  and  Case  4  are  similar  to  each  other  with  more  lateral  strain  in 
Case  4.  This  is  consistent  with  the  material  being  harder  in  the  TD  than  in  the  RD  as  shown 
in  the  uniaxial  tests  (see  Figure  40  which  show  that  for  strain  levels  below  15%,  the  material 
is  stronger  in  the  TD  direction  for  both  tensile  and  compressive  loading).  For  comparison 
purposes,  FE  simulations  using  the  isotropic  von  Mises  yield  criterion  in  conjunction  with  an 
isotropic  hardening  law  based  on  the  RD  tension  data  were  conducted  and  compared  to  the 
results  obtained  using  the  anisotropic  model  (see  Figure  41).  Note  that  in  Case  1  and  Case 
3,  the  hard  direction  (TT)  is  the  width  direction  and  the  isotropic  simulation  shows  more 
deformation  than  it  is  obtained  using  the  anisotropic  model. For  Case  2  and  Case  4  there  is  less 
deformation  in  the  height  of  the  beam  and  the  width  shows  deformation  similar  to  the  isotropic 
simulation. Again  this  is  because  for  these  loading  scenarios,  the  softer  direction  is  aligned  with 
the  width.  For  each  loading  scenario,  comparisons  between  FE  axial  strain  contours  obtained 
with  the  proposed  anisotropic  model  and  experimental  axial  strain  contours  obtained  by  digital 
image  correlation  are  presented.  The  beam  orientation  for  Case  1  loading  is  shown  in  Figure 
42  for  reference  while  Figure  43  shows  the  comparison  of  the  respective  strain  profiles.  A  very 
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Case  1 


Case  2 


Figure  39:  Simulated  final  cross  sections  for  the  four  beam  configurations 


good  agreement  is  obtained.  Note  that  the  data  from  the  experiment  does  not  cover  the  entire 
profile  area  due  to  the  DIC  [[42,  43]]  techniques  used. 

A  more  qualitative  comparison  is  presented  in  Figure  44,  which  shows  a  plot  of  the  axial  strain 
versus  the  height  of  the  beam  at  the  center  of  the  beam.  A  very  good  agreement  between  the 
experiment  and  simulation  is  observed.  In  particular, the  upward  shift  of  the  neutral  axis  of 
the  beam  is  very  well  captured.  To  provide  a  final  validation  of  the  model,  the  beams  were 
sectioned  at  the  midpoint  and  a  photograph  of  the  cross  section  was  compared  to  the  simulation 
results. The  comparison  between  theory  and  data  for  Case  1  is  shown  in  Figure  45.  Note  that 
there  is  very  little  deformation  perpendicular  to  the  loading  direction  because  this  is  the  harder 
to  deform  direction. 

The  beam  orientation  for  Case  2  is  shown  in  Figure  46  for  reference,  while  Figure  47  shows 
the  comparison  of  the  respective  strain  profiles.  Again,  very  good  agreement  is  shown  between 
experimental  data  and  the  simulation  results  obtained  using  the  proposed  model.  Note  that 
the  model  predicts  somewhat  less  strain  through  the  thickness  in  the  loading  direction. 
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Figure  40:  Comparison  between  stress-strain  response  in  tension  and  compression  and  experi¬ 
mental  data  for  RD  and  TD  directions,  respectively. 


Case  2 


Case  4 


Figure  41:  Comparison  of  the  final  cross  sections  of  the  beam  specimens  as  predicted  by  the 
proposed  model  and  by  an  isotropic  von  Mises  model  (red  mesh) 
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Figure  42:  Case  1:  Long  axis  in  RD  direction,  loading  in  the  TD  direction. 


x  (mm) 


Figure  43:  Case  1  loading:  comparison  of  FE  axial  strain  contours  (ex)  obtained  with  the 
proposed  model  against  experimental  data:  £=RD,  y=TD. 
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Figure  44:  Plate  1,  Case  1:  Axial  strains  ( ex )  versus  height  at  centerline:  :i;=RD,  j/=TD. 
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Figure  45:  Case  1:  Comparison  of  the  photographed  experimental  cross  sections  and  FE  simu¬ 
lation  using  the  proposed  model(symbols):  y= TD,  z= TT. 


Figure  46:  Case  2  loading:  Long  axis  in  the  RD,  loading  in  TT 
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Figure  47:  Case  2:  Comparison  of  the  photographed  experimental  cross  sections  and  FE  simu¬ 
lation  using  the  proposed  model(symbols):  :/;=RD.  z=TT 


Figure  48:  Case  2:  Axial  strains  (ex)  versus  height  at  centerline:  :/;=RD.  z=TT 
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A  plot  of  the  axial  strain  versus  the  height  of  the  beam  at  the  center  of  the  beam  is  shown 
in  Figure  48.  Note  the  very  good  agreement  between  the  experiment  and  simulation  and  the 
clear  upward  shift  of  the  neutral  axis  of  the  beam  due  to  the  tension-compression  asymmetry 
of  the  material,  which  is  very  well  described  by  the  model. 

The  comparison  of  final  cross  sections  of  the  beam  in  Case  2  loading  is  shown  in  Figure 
49.  Again,  there  is  a  very  good  agreement.  Note  that  the  model  predicts  with  accuracy  that 
there  is  more  deformation  perpendicular  to  the  loading  direction  because  in  Case  2,  the  loading 
direction  is  the  softer  transverse  direction.  The  beam  orientation  for  Case  3  is  shown  in  Figure 


Figure  49:  Case  2:  Comparison  of  the  photographed  experimental  cross  sections  and  FE  simu¬ 
lation  using  the  proposed  model  (symbols):  ?/=TD,  z=TT 


50  for  reference  while  Figure  51  shows  the  comparison  between  experimental  and  simulated 
strain  profiles.  Again,  very  good  agreement  is  obtained. 

A  plot  of  the  axial  strain  versus  the  height  of  the  beam  at  the  center  of  the  beam  is  shown 
in  Figure  52.  Note  the  very  good  agreement  between  the  experiments  and  simulation  results. 
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Figure  50:  Case  3:  Long  axis  in  the  TD  direction,  loading  in  the  RD  direction 


The  upward  shift  of  the  neutral  axis  is  very  well  described  by  the  model.  The  comparison  of 
the  photographed  and  simulated  final  cross  sections  is  shown  in  Figure  53.  Note  the  excellent 
agreement. 
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Figure  51:  Case  3:  Comparison  of  the  FE  axial  strain  contours  (ey) obtained  with  the  proposed 
model  against  experimental  data:  x=RD,  y—TD 
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Figure  52:  Case  3:  Axial  strains  (ey)  versus  height  at  centerline:  :r=RD,  y=TD 
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Figure  53:  Case  3:  Comparison  of  the  photographed  experimental  cross  sections  and  FE  simu¬ 
lation  using  the  proposed  model(symbols):  .x=RD,  z=TT 
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The  beam  orientation  for  Case  4  is  shown  in  Figure  54  for  reference  while  Figure  55  shows 
the  comparison  of  the  simulated  and  experimental  axial  strain  contours  for  Case  4  .  Again, 
very  good  agreement  is  shown. 


z  =  TT 


x  =  RD 


y  =  TD 


Figure  54:  Case  4:  Long  axis  in  the  TD  direction,  loading  in  the  TT  direction 
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Figure  55:  Case  4:  Comparison  of  experimental  and  simulated  axial  strain  contours  (ey)  : 
y=TD,  z= TT 


The  axial  strain  at  the  center  of  the  beam  along  the  height  of  the  beam  is  shown  in  Figure  56. 
Note  the  very  good  agreement  between  the  experimental  and  simulation  results,  in  particular 
the  upward  shift  of  the  neutral  axis  of  the  beam,  which  is  due  to  the  tension-compression 
symmetry  of  the  material.  A  very  good  agreement  between  experimental  and  simulated  final 
cross  sections  for  Case  4  is  also  obtained  (see  Figure  57). 
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Figure  56:  Case  4:  Axial  strains  (ey)  versus  height  at  centerline:  y=TD,  z=TT 


Figure  57:  Plate  1,  Case  4:  Comparison  of  the  photographed  experimental  and  simulated 
(symbols) cross  sections  :  x=RD,  z=TT 
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4  Characterization  of  the  high-strain  rate  behavior  of 
high-purity  o-titanium 


High  rate  loading  tests  were  carried  out  using  a  split  Kolsky-Hopkinson  pressure  bar  at  Eglin  Air 
Force  Base,  Florida.  High  rate  compression  specimens  (Figures  58  A)  )  were  simple  cylinders 
of  0.2  x  0.2  inches  thus  smaller  than  the  specimens  used  for  the  quasi-static  compression  tests. 
High-rate  tensile  specimens  were  cylindrical  dog  bones  as  shown  in  Figure  58  (B)  and  were 
marked  with  an  arrow  indicating  the  top  of  the  plate  and  therefore  the  through  thickness 
direction.  The  plate  thickness  did  not  allow  enough  material  to  obtain  high-rate  through 
thickness  tension  specimens. 


A)  B) 


Figure  58:  High  rate  test  specimens:  A)  for  compressive  loading;  B)  for  tensile  loading. 

All  tensile  tests  were  carried  out  to  failure.  A  typical  cross  section  of  the  failed  specimen 
is  shown  in  Figure  59. Note  the  cup-cone  failure  and  the  elliptical  shape  of  the  final  cross- 
section,  with  the  harder,  through-thickness  direction  in  the  direction  of  the  major  axis  of 
the  specimen.  Both  compressive  and  tensile  tests  were  carried  out  on  specimens  cut  along 
the  rolling,  transverse,  and  through-thickness  directions, respectively  at  nominal  strain  rates 
ranging  from  400  to  600/s  (see  Fig.  60  and  Tables  5  and  6  ).  Table  7  gives  the  yield  values 
for  several  levels  of  strain  as  well  as  the  anisotropy  ratio  (defined  as  the  ratio  of  highest  to 
lowest  yield  at  a  given  strain  level)  for  both  the  high  loading  rate  data  and  the  quasi-static 
data. The  test  results  show  that  the  material  response  at  high-strain  rate  is  also  orthotropic. 
The  transverse  direction  stress-strain  curve  remains  below  the  through-thickness  stress-strain 
curve  throughout,  as  in  the  quasi-static  case.  As  in  the  case  of  quasi-static  loading,  the  plate 
material  is  initially  harder  in  the  through-thickness  direction  but  after  15%  strain,  the  flow 
stress  in  the  rolling  direction  is  the  largest. 
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Side  View 


Figure  59:  Post-test  specimen  with  long  axis  along  in-plane  direction,  subjected  to  high-strain 
rate  tensile  loading  showing  cup-cone  fracture  and  elliptical  final  cross-section. 


Table  5:  Strain  rates  achieved  in  tensile  Split-Hopkinson  Bar  (SHPB)  tests 


Test  Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Strain  Rate  (sec-1) 

522 

517 

643 

665 

652 

627 

662 

653 

555 

543 

Test  Number 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

Strain  Rate  (sec-1) 

585 

534 

547 

563 

573 

528 

552 

557 

529 

516 

Test  Number 

21 

22 

23 

Strain  Rate  (sec-1) 

517 

616 

531 

Note:  AVG  =  567  and  Standard  Dev  =  53 


_ Table  6:  Strain  rates  achieved  for  tensile  SHPB  tests _ 

Test  Number  123456789  10 

Strain  Rate  (sec”1)  419  429  318  369  396  466  457  408  417  403 

Note:  AVG  =  407  and  Standard  Dev  =  42 
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Stress  (MPa) 


Figure  60:  Experimental  compression  results  showing  that  the  high- strain  rate  response  at 
400/s  is  orthotropic. 
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Table  7:  Quasi-static  and  high  rate  compressive  yield  values  in  the  RD,  TD,  and  TT  directions, 
respectively  and  anisotropy  ratios  for  the  material  under  high-rate  loading 


High  Rate 

Quasi-static 

Strain 

RD 

TD 

TT 

Max/ Min 

RD 

TD 

TT 

Max/ Min 

0.05 

411 

398 

294 

1.398 

225 

276 

330 

1.467 

0.1 

498 

467 

399 

1.248 

271 

307 

361 

1.332 

0.15 

537 

498 

492 

1.091 

323 

331 

389 

1.204 

0.2 

548 

521 

573 

1.099 

379 

352 

419 

1.190 

0.25 

589 

548 

626 

1.142 

426 

388 

448 

1.155 

0.3 

620 

575 

654 

1.137 

456 

418 

477 

1.141 

0.35 

NA 

NA 

NA 

NA 

493 

452 

505 

1.117 

0.4 

NA 

NA 

NA 

NA 

522 

483 

537 

1.112 

Note  also,  the  strongly  non-linear  strain  hardening  behavior  of  the  material  in  the  RD  com¬ 
pression,  with  an  abrupt  change  in  the  slope  of  the  stress-strain  curve  at  15%  strain,  which  is 
indicative  of  significant  twinning.  The  change  in  hardening  rate  is  even  more  dramatic  at  high 
rate  than  at  quasi- static  loading  rate.  This  is  not  complectly  unexpected  as  it  has  been  ob¬ 
served  by  [1]  that  titanium  twins  more  readily  as  loading  rate  increases.  In  general,  irrespective 
of  the  loading  orientation,  the  yield  stress  increases  with  increasing  loading  rate. 

The  tests  have  also  revealed  that  the  tension-compression  asymmetry  in  the  high-strain 
rate  regime  is  highly  directional.  Figure  61  A)  shows  a  comparison  between  the  high  -rate 
results  for  uniaxial  loading  in  the  TD  direction  in  tension  and  compression,  respectively.  The 
initial  yield  points  are  very  close  but  as  more  deformation  occurs  the  strength  in  compression 
becomes  somewhat  larger.  This  may  indicate  that  some  twinning  is  occurring  in  the  compression 
loading  but  this  has  not  been  confirmed  by  post  test  metallography.  Comparisons  of  high  rate 
compression  data  with  quasi-static  data  for  the  TD  are  shown  in  Figure  61  B).  A  clear  increase 
in  strength  is  observed  but  the  hardening  rate  remains  nearly  unchanged. 

The  high  rate  results  for  compressive  uniaxial  loading  in  the  TT  direction  compared  to  data 
gathered  at  quasi- static  loading  rates  is  shown  in  Figure  62.  Again,  a  clear  increase  in  strength 
with  loading  rate  is  observed  with  very  little  change  in  hardening  rate.  Due  to  geometry 
constraints,  no  through  thickness  tension  data  could  be  obtained  at  high  rates  of  loading. 

Results  at  high-  rate  uniaxial  loading  in  the  RD  are  shown  in  Figure  63  A).  As  for  loading 
in  TD  direction,  the  initial  yield  levels  are  similar  but  the  difference  increases  with  additional 
strain.  This  difference  increases  at  a  much  higher  rate  than  for  the  TD  data.  There  is  also  a 
clear  change  in  hardening,  which  indicates  that  significant  twinning  is  occurring.  Figure  63  B) 
shows  the  comparison  with  quasi-static  data  where  post  test  metallography  showed  a  significant 
amount  of  twinning.  The  data  indicates  that  even  higher  levels  of  twining  may  be  occurring  at 
the  higher  loading  rates. 
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A  B 

Figure  61:  A)  Experimental  high  rate  data  for  uniaxial  loading  in  the  TD  direction  in  tension 
and  compression,  respectively,  and  B)  Effect  of  the  strain-rate  in  the  response:  comparison  of 
experimental  high  rate  data  to  experimental  quasi-static  data 


Figure  62:  Comparison  of  compressive  high  rate  to  quasi-static  data  in  the  TT  direction 
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Figure  63:  A)  Comparison  between  experimental  high  rate  response  in  uniaxial  loading  along 
the  RD  direction  for  tension  and  compression,  respectively;  B)  Comparison  of  high  rate  data 
to  quasi-static  data  for  loading  along  the  RD  direction 
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5  Discussion  and  Conclusions 


As  already  stated  ,  a  primary  goal  of  this  research  was  to  advance  the  current  state-of-the-art 
by  developing  user-friendly,  micro-structurally  based  and  numerically  robust  macroscopic  con¬ 
stitutive  models  that  can  capture  with  accuracy  the  particularities  of  the  plastic  response  of 
hexagonal  metals,  in  particular  high  purity  titanium.  It  has  been  demonstrated  that  this  goal 
has  been  met  to  a  large  degree.  A  comprehensive  set  of  experiments  were  conducted  in  or¬ 
der  to  understand  the  response  of  high-purity  titanium  at  room  temperature  and  consequently 
model  the  observed  response  at  different  strain  rates.  A  series  of  monotonic  quasi-static  uni¬ 
axial  tension  and  compression  tests  were  conducted  to  quantify  the  plastic  anisotropy  and  the 
tension-compression  asymmetry  of  this  material.  In  addition,  four-point  bending  tests  were 
conducted.  To  characterize  the  material’s  strain  rate  sensitivity,  Split  Hopkinson  Pressure  Bar 
tests  at  strain  rates  of  400  to  600  sec-1  were  performed  along  the  axes  of  symmetry  of  the 
material.  To  gain  insight  into  the  deformation  mechanisms,  OIM  and  texture  evaluations  were 
also  conducted  in  conjunction  with  the  mechanical  tests.  The  experimental  activities  were 
complemented  with  multi-scale  model  development  in  the  framework  of  plasticity.  The  ma¬ 
terial  parameter  identification  procedure  is  also  outlined.  The  algorithmic  aspects  related  to 
the  implementation  of  the  models  developed  in  the  Finite-Element  code  EPIC  are  presented 
in  detail.  The  ability  of  the  proposed  model  to  capture  the  main  features  of  the  observed 
behavior  is  examined  by  comparing  the  experimental  data  with  simulation  results  in  terms  of 
stress-strain  response  in  tension/compression  and  bending  for  different  orientations  with  re¬ 
spect  to  the  axes  of  symmetry  of  the  material.  Comparison  between  proposed  theory  and  data 
are  extremely  good.  This  integrated  experimental-theoretical  effort  appears  to  have  resulted  in 
the  description,  for  the  first  time,  of  the  anisotropic  stress-response  of  high-purity  o-titanium 
at  room  temperature.  Further  research  is  needed  to  explore  other  loading  environments  but 
this  work  has  shown  that  the  proposed  model  can  describe  with  great  accuracy  the  anisotropy 
induced  by  texture  evolution  during  monotonic  loading.  The  ability  to  incorporate  data  from 
other  loading  conditions  is  already  in  place.  Although,  this  research  was  concerned  primarily 
with  high  purity  titanium,  it  is  clear  that  the  proposed  model  and  implementation  approach  is 
applicable  to  other  hexagonal  materials. 
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